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Abstract 

Statistical method is applied to Neveu and Schwartz model for obtain- 
ing the average characteristics of the heavy resonances: spin distribu- 
tion, decay widths etc. The properties of the dual model spectrum of 
states constructed of both commuting and anticommuting operators 
are considered. In this case the spin characteristics coincide with the 
results of the statistical bootstrap model. 
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1 Introduction 



In this paper we investigate the average characteristics of the Neveu and 
Schwarz dual model from the statistical view-point [1 — 4]. The exponential 
rise of the density of hadron states at high energies in the dual resonance 
model (DRM) creates the preconditions for the efficiency of the statistical 
approach. Besides, there may exist an essential connection [3] between the 
properties of the resonance spectrum in the DRM and in the statistical boot- 
strap model [5, 6]. 

For some observable A the average in the resonances with fixed mass M 
(i.e. microcanonical average) a(M 2 ) = a(0) + a'M 2 = N is defined in dual 
resonance model to be [3] 

(A) ^N imR) ±§dzZ-»-^{ Z »A) 

{ )N £ N ( R \ R ) ±fdzz-*-i8j>z* ' 1 } 

where H is the Hamiltonian operator in DRM, and the integration contour 
in the complex z-plane envelopes the point z — 0. 

In the paper [4] a new type of averaging has been proposed 



E \(0\V(p)\R)\ 2 (R\A\R) 

((A)) » = 



E \(0\V(p)\R)\ 
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^fdxfdyx- N - 1 y- N - 1 (0\V(p)x H Ay H V+(p)\0) 

±fdxx-N-i(0\V(p)x*V+(p)\0) ' (2) 

where V(p) is the vertex operator in DRM. Every resonance state is contained 
in the average (2) with the weight proportional to the square of the coupling 
constant with the basic states of the model (conditionally 7r-mesons) g n7r R = 
(0\V(p)\R). This provides an effective accounting of the dynamical features 
of the reactions and allows us to call (2) the dynamical average. 

In the papers [3, 4] the averages such as (1) and (2) were considered in 
the framework of the Veneziano dual model. It is clear that the statistical 
approach can be applied to the other possible formulations of the dual model. 
We shall consider here the dual model of Neveu and Schwarz which possesses 
some important theoretical advantages comparing with the usual Veneziano 
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model. In this model the resonance states are generated by the infinite 
sequence of operators of the usual Veneziano model 

K,<n = -<r<w (3) 

(n,ri = 1,2,...) 
and the new infinite set of anticommuting operators 

{b^Kt} = ~9^ mm , (4) 

(m,m' = 1/2,3/2,...) 

satisfying conditions [a£, &™ + ] = 0. The role of Hamiltonian is played by the 
operator 

oo oo 

H = H a + H b = -J2 na\a n - mb mb m (5) 

n=l m=l/2 

and the vertex operator is defined as 



V(p) = V a (p)V b (p) = exp [ - ^Et] ex P [^Et1^(^ + bm). 

(6) 

The eigenstates of the operator H are divided into two groups with different 
quantities of G-parity [7] (G = (-l)S- fe ™ bm ); 

1) the states with eigenvalues N — |, |, ... composed of the odd number 
of "6" -operators and an arbitrary number of "a" -operators, G = — 1; 

2) the states with eigenvalues N — 1,2, ... composed of the even number 
of "b" -operators and an arbitrary number of "a" -operators, G — 1 (the 
state with N = is spurious [7]). 

In the particle interactions the value of G-parity is conserved. The 7r- 
meson is associated with the lowest physical state of the model with the 
mass ml = —\. 

The low energy part of the spectrum in the Neveu-Schwarz model is 
similar to the experimental data. 
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2 Statistical Averages 

We shall start with considering the degeneracy of the resonance level N, i.e. 
the number d(N) of the solutions of 

H\R i )=N\R i ),(i = ^,---d(N)). (7) 
It is obvious that the degeneracy d(N) can be presented in the following form 



d<K N) = _L fdzz-™-^? (z 2H ) = — Idzz'™' 1 x 

OO I 

x rf(l -z 2n )-\l + z 2n - 1 ) A = — i dzz-^-^iOMefiOWUz 1 / 2 , (8) 
„ i 2m J 

where r = — and 62, 4 are the Jacobi ^-functions [8]. The behaviour of 
(8) as iV goes to infinity is determined by the quantity of the integrand at 
z ~ 1, which is given by means of the imaginary Jacobi transformation [8]. 
The integral (8) is calculated by employing the saddle-point technique, and 

T=-VN (9) 

7T 

determines the saddle point Z = e~^r . We get 

d(N) ~ ^N~ 7/4 exp[2nVN}. (10) 

Taking into consideration that N ~ a'M 2 , we get the expression for the 
hadron state density p(M) at M — > 00 



p(M) ~ const • M 5/2 exp 



M 



T, 



(11) 



where To = 2 J^j plays the role of Hagedorn limiting temperature [5]. For 
a! ~ 1 GeV~ 2 we get To ~ 160 MeV, and this result remarkably coincides 
with the prediction of the thermodynamic model [5] based on the comparison 
with the experiment. 

Let us calculate now the average number of the oscillator excitations. 
According to (1), we have 
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{~a + n a n ) N = — f dz z-™-^ ( - a+ n a n z™) 



d{N)2m J l-z 2n 



dzz-™- 1 ^-^ n a - * 2 T 4 a + ^t- (12) 

fc=l 



At fixed n and A" — > oo we get 

4 

(-a+a n )jv ~ e „ /r _ 1 , (13) 

where T is given by (9). 

Analogous calculation for the number of "6" -excitations gives 

(-btb m ) N * (14) 

Expressions (13) and (14) are the equilibrium distributions of Bose-Einstein 
and Fermi-Dirac, respectively. Thus, we have come to the canonical distri- 
bution with the temperature T given by (9). The use of the saddle point 
technique while calculating the integrals such as (1) enables us to fulfil the 
normal in statistical mechanics transition from the microcanonical to the 
canonical (Gibbs) ensemble. The canonical averages of the physical values 
are defined as 

{A)T= Spe-^ • (15) 

Let us continue considering the average characteristics of the resonance spec- 
trum using (15). The total mean number of the oscillator excitations of "a" 
and "6" type are equal, respectively, to 

(- E «Wt = E^rT^ 4T ln T > (i6) 

n=l,2,... n e 1 

(- E ^UT = E^TT" 4Tln2 - ^ 

m=l/2,3/2,... rn ~r 1 

The canonical average of energy has to coincide with the energy A" of the 
microcanonical ensemble and is given by the expression: 
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(H) T = Wt+ Wt = £ ^ = (18) 

In (16), (17), (18) we substitute the summation for the integration which is 
reasonable for large T. 

Let us proceed to the question of spin distribution in the Neveu and 
Schwarz dual model. In the case of Veneziano model this task has been 
solved in the papers [2, 9]. The operator of the z-component of spin is [10] 

00 oo 

l, = li + l\ = iY,K + < - < + <) + * E (C^ - TO- (19) 

n=\ m=l/2 

Choosing the basis of states to be diagonal in L z [2] , we find the characteristic 
function of the operator L z according to (15) 

1 — n IT\ 2 1 , r> n-1/2 2n-l 

l — e n/1 ) z i r 1 + 2 cos (pe t +e t -, 



(e itfiL *) T = TT [ - - - 

„=i 1 - 2 cos ¥>e- n / T + e~ 2n / T 



(1 + e — ^) 2 



n=l 

= 9^ 1 (v\t)9 2 (0\t)9 3 (v\t)9 4 (0\t) sin™, (20) 

where v = r = and 9\, 62, #3, #4 are the Jacobi ^-functions. For 
T — > 00 we apply the imaginary Jacobi transformation for the evaluation of 
(20), and get 

, ■ r . 2nT sin ^ . , 

^ = (21) 

By differentiating (21) we find the average value 

H 2 tt 2 T 2 N 

The distribution in Z 2 is the following 
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The integral in (23) is calculated by contour integration enclosing the con- 
tour in the lower half-plane and explicitly performing the sum over the pole 
residues. Using (23) one can easily obtain the spin distribution of the reso- 
nances 

s d , , , 1 (21 + 1) sinh (4f) , s 

MI) . _ (a + d w mq\^ = w : js^r - (24) 

which can give further information about the spin structure of the resonance 
spectrum. Let us write out the spin characteristics of the resonances con- 
structed of "6" -operator only. The characteristic function of the operator L b z 
is 



(n-l/2) (2w-l) 

( >T "ii (l + e-^) 2 -^ 3 (0|r)T-^ e • ^ 

From (25) we get 

^T=~(^ Li ) T \^^T=^. (26) 

We see that for N — > oo the "6" -oscillator contribution into spin becomes 
asymptotically negligible, whereas their contribution into the energy makes 
up y according to (18). The distribution in Z 2 for the "6" oscillators has the 
Gaussian form 

*<'■> = s£^ e -* T *7£r e ~*- (27) 
3 Dynamical Averages 

The dynamical averaging (2) in Neveu-Schwarz model is performed over the 
states with integer values of N (it is necessary for the G-parity conservation). 
The characteristic function of the operator L z is given by 

£-Jdxx- N -'(i;\V(p)x H e^V+(p)\n) 
* ))N ±fdxx-»-HK\V( P )x*V+(p)\n) ' 1 ] 
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where \ir) = kby 2 \0) is pion state [7]. Using kinematical relations in the 

c.m.s. (p + k = 0) and the mass-shell conditions (p 2 = k 2 = m 2 = — |), we 
obtain 



(n\V(p)x H e lipL *V + (p)\ir) = [(1 - x )-W(pI+p 2 v )(i-™ <p)} x 
x[\y^ + ^(1 -x) + (pl+p 2 y )(Nx - T ^-)(l - cos*,)]. (29) 

From (28) and (29) we find 



((L 2 z ))n = -£-2^ Lz ))n\^o * KpI + PDIHN -1) + C + 0(1)], (30) 

where C — 0, 577 ... is the Euler constant. 

Choosing the momenta in the c.m.s. along the x-axis, for (30) we get 

/V 

((L 2 z )) N ^-\nN. (31) 
Also, it is easy to find the contribution of "6" -oscillators 

(((L b z ) 2 }) N ^4(p 2 x +p 2 y )^N. (32) 
The average numbers of the excitations are 

((-a+a n ))iv = -(1 - ^);forn < N, (33) 
{{-b + m b m )) N = ^-M\<m<N. (34) 



4 Resonance Decay Widths 

Let us take up now a heavy resonance decay within the framework of the 
Neveu and Schwarz model. According to [3], the average decay width of the 
resonance from iVi-level into the resonance from N2- level and 7r-meson is 
given by 
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<r> 



1 



d(N 1 )(2TTi) 2 

The trace in (35) has the following form 



jdx jdyx- 2N '- l y- 2N *- l $v{x 2H V(p)y 2H V + (p)). (35) 



Sp (x 2H V{p)y 2H V + {p) 
(1 - (xy) 2n ^ 2 



[1 - X 2 )' 1 H(l - (xy) 2n )- 4 x 



71=1 



X 



1 -x 2 (xy) 2n )(] 



oo x 2k-l + y 2k-l 



^^][in (1+( ^-r £ ^^]. (36) 



To calculate the integral (35), we introduce the variables uj = xy and x. The 
main contribution to the integral (35) with large N 2 , Ni — N 2 is given by 
the region uj ~ l,x ~ 1. Expressing the integrand in terms of the Jacobi 
^-functions, we obtain for uj ~ 1 



n (T 

71=1 I 1 



{l-uj 2n f 



7r sinh(lna: 



r ln 2 1, 



x 2 uj 2n )(l - ±uj 2n ) lnwsin(^f) eXPL lnw J ' 



oo x 2k-l 



E 

k=i 



+ (*) 



2fc-l 



7T 



1 + CU 



2fc-l 



2 In uj sin (^V 

v In a; / 



(37) 



(38) 



The integration over uj is performed by employing the saddle-point technique 
with the result 



d(N 2 ) 1 



d(N ± ) 2ni 

7T 



jdxx~ 2N ^ 2N *-\l - x^fix,^); 
2 (— sinh In a;) 



x^i 21nc<j sin 



,2 7rlri3: 
In wo 



(39) 



where ujq = exp 



, j is the saddle point. The remaining integral is easily 
calculated at N\ — N 2 > > V^2 which corresponds to the large energy of the 



emitted 7r-meson. In this case 



Inx 



/N2 



In uiq 



Ni-N 2 



0, and so we get 



<r> 



d(N 2 ) Ni - N 2 
d(Ni) 2 



(40) 
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Expression (40) may be written in terms of the decaying resonance mass Mi 
and the 7r-meson energy E 



<r> ^ (i 



2E 



Y 7/A EM 1 exp [ 



(41) 



Mi 



We see that the meson energy spectrum has the Boltzman distribution form 



with the effective temperature T e ff = -^-(1 + y 1 — j^). Such a dependence 

coincides exactly with the results of the statistical bootstrap model [12] (note 
that we do not take into account the whole chain of the successive decays, 
the consideration of which changes the value of T e ff [12]). 

A considerable difference of (41) from the analogous result obtained in 
the Veneziano model with a(0) ^ 1 [3] consists in the linear growth of the 
decay width with Mi increasing (such a behaviour would take place also in 
the Veneziano model with a(0) = 1). 

5 Discussion 

In conclusion we would like to discuss some interesting, as it seems to us, 
aspects of our consideration: 

1. The average resonance spectrum characteristics in the Neveu-Schwarz 
model for iV — > oo are similar, as a whole, to those obtained in the Veneziano 
model [2, 3, 4]. This allows us to hope that statistical approach grasps the 
main features of the dual dynamics independent of the choice of the specific 
model. 

2. The average decay width has the same dependence on the 7r-meson 
energy as in the statistical bootstrap model and grows linearly as M — > oo. 
Concerning the total width dependence on the hadron mass, the comparison 
is meaningless, since no decay widths but the relative probabilities are con- 
sidered in the statistical bootstrap and the total probability is normalized to 
be equal unit. 

3. In the Neveu-Schwarz model there exists an infinite sequence of the 
gauge relations which allows one to remove the unphysical states. In order to 
obtain the number of physical states, we must add the factor (1 — z) n^Li (1 _ 
z 2n )(l + 2 2n ~ 1 )~l to the integrand (8) (the same factor as for calculating the 
planar loop) [13, 11]. The number of physical states is given by [13] 
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d D Phys {N) = d D -\N)-d D -\N- l -) ) (42) 

where D is the dimension of the oscillator operators (in our case D — 4). For 
large M the density Pp^ ys equals 

pPhys(M) ~ constM" 3 exp[7Tv / 3o 7 M]. (43) 

The removal of unphysical states can be easily taken into account in our 
calculations. For example, one must substitute (42) for d(Ni) and d(N 2 ) in 
the expression (40), which results in the alteration of the quantity of To in 
(41). So, the partial widths of the decay with 7r-meson emission, which have 
the form (T) ~ , do not essentially change as the removal of unphysical 
states is performed. However, the situation may be quite different while 
cosidering the general case of decay Ri — > R 2 + R3 where one expects the 
results such as (r) ~ ^^^j^- Here the decay characteristic will be crucially 
dependent of the value of power of M in the expressions (11), (43). Note that 
the case of the power value equal to —3 (regarding the physical states only) 
coincides with the results of the investigations [14, 15] and makes the decay 
into one heavy and one light (tt- meson) particles dominant, while the case of 
the value equal to —5/2 (all states are taken into account) corresponds to the 
primary formulation of Hagedorn [5] which predicts the decay into particles 
of approximately equal masses. 

4. Our calculations allow us to obtain the average characteristics of the 
resonance spectrum constructed of "6" operators only. 

The question of building the dual theory with the same hadron state 
density as in the statistical bootstrap has been considered in the recent paper 
[16]. The author concludes that the construction of resonance states with the 
only Fermi operators represents one of the most preferable possibilities. In 
this case using calculations (26) and (27) we find that also the average spin 
and the spin distribution coincide with that of the statistical bootstrap theory 
with spin [17]. 

We thank B. V. Struminsky and V. A. Miransky for useful discussions. 
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